GLOBAL WEAK SOLUTIONS TO THE INHOMOGENEOUS 
NAVIER-STOKES-VLASOV EQUATIONS 
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Abstract. A fluid-particle system of the incompressible inhomogeneous Navier-Stokes 
equations and Vlasov equation in the three dimensional space is considered in this paper. 
The coupling arises from the drag force in the fluid equations and the acceleration in 
the Vlasov equation. An initial-boundary value problem is studied in a bounded do- 
main with large data. The existence of global weak solutions is established through an 
approximation scheme, energy estimates, and weak convergence. 



1. Introduction 

We are concerned with the Navier-Stokes-Vlasov equations for particles dispersed in a 
density-dependent incompressible viscous fluid: 

pt + div(pu) = 0, (1.1) 

(pu) t + div(/9U <g> u) + Vp - pAu = - / m p Ffdv, (1.2) 

Jr 3 

divu = 0, (1.3) 

/t + v-V a; / + div v (F/) = ) (1.4) 

for (x, v, t) in O x R 3 x (0, oo), where C M 3 is a bounded domain, p is the density of the 
fluid, u is the velocity of the fluid, p is the pressure, p is the kinematic viscosity of the fluid. 
The density distribution function f(t,x,v) of particles depends on the time t € [0, T], the 
physical position i £ !! and the velocity of particle v 6 ]R 3 . In (jl.2p . m p is the mass of 
the particle and F is the drag force. The interaction of the fluid and particles is through 
the drag force exerted by the fluid onto the particles. The drag force F typically depends 
on the relative velocity u — v and on the density of fluid p (e.g. |18j). such as 

F = F p(u-v), (1.5) 

where Fq is a positive constant. Without loss of generality we take p = Fq = m p = 1 
throughout the paper. The fluid-particle system f)l. .5[) arises in many applications 
such as sprays, aerosols, and more general two phase flows where one phase (disperse) can 
be considered as a suspension of particles onto the other one (dense) regarded as a fluid. 
This system (|l.ip - ([1.5p or its variants have been used in sedimentation of solid grain by 
external forces, for fuel-droplets in combustion theory (such as in the study of engines), 
chemical engineering, bio-sprays in medicine, waste water treatment, and pollutants in the 
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air. We refer the readers to [TJ 0J [9j [101 El El [201 [22] for more physical background, 
applications and discussions of the fluid-particle systems. 

The aim of this paper is to establish the global existence of weak solutions to the 
initial-boundary value problem for the system (|l.ip - (|1.5p subject to the following initial 
data: 

p|t=o = Po(x) > 0, x G fi, (1.6) 
(pu)|t =0 = mo(i), x G Q, (1.7) 
/|t=o = /o(s,v), xefi, v G M. N , (1.8) 
and the following boundary conditions: 
u(i, x) = on dil, 

(19) 

f(t, x, v) = f(t, x, v*) for x G <9ri, v • < 0, V ' 7 

where 

v* = v — 2(v • z/(x))z/(x) 

is the specular velocity, and fix) is the outward normal vector to ft. The same as in |14j . 
we assume that the particles are reflected by the boundary following the specular reflection 
laws. When the drag force is assumed independent of density in (|1.5h . the hydrodynamic 
limits and the global existence of weak solutions to the Navier-Stokes and Vlasov-Fokker- 
Planck equations were studied in [121 El El E] • When the drag force depends on the 
density as in (jl.5p . a relaxation of the kinetic regime toward a hydrodynamic regime 
with velocity u on the vacuum {p = 0} can not be excepted. It is hard to establish 
a priori lower bounds on the density from the mathematics view point. The objective 
of this paper is to establish the existence of global weak solutions to the initial-value 
problem (|l.ip - (|1.9p with large data in certain functional spaces. In general, the analysis 
of fluid-particle system is challenging since the density distribution function / depends 
on more variables than the fluid density p and velocity u. The existence of global weak 
solutions to the Stokes- Vlasov equations in a bounded domain was studied in [13]. In [2] 
the convection term was included and the incompressible Navier-Stokes- Vlasov equations 
were considered in periodic domains. A similar system with thermal diffusion acting 
on the particles, that is, the incompressible Navier-Stokes equations coupled with the 
Vlasov- Fokker-Planck equations, has been studied in [11], where the authors established 
the existence of classical solutions with small data. Recently, the existences of global 
solutions to the incompressible Navier-Stokes- Vlasov equations in a bounded domain or in 
the whole space were obtained in [231 124] . Also there are a lot of works on hydrodynamic 
limits, we refer the reader to [U [121 [131 [18] (and the references therein) where some 
scaling and convergence methods such as the compactness and relative entropy method 
were applied to investigate the hydrodynamic limits. A key idea in [121 [13] is to control 
the dissipation rate of a certain free energy associated with the whole space. The global 
existence of weak solutions to the compressible Navier-Stokes equations coupled to Vlasov- 
Fokker-Planck equations was established in [17]. The coupled system (ll.l| )- (jl.5p has extra 
difficulties due to the appearance of density in the interactions and in the Vlasov equation 
as well as the lack of diffusion in the Vlasov equation. We note that the local classical 
solution to the Euler- Vlasov equations was obtained in [TJ when the drag force F is assumed 
to be in the form of (II. 5p . 
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When / does not appear, the system (|l.ip - (|1.4p is reduced to the density-dependent 
Navier-Stokes equations. We refer readers to Lions [15] for the compactness, the exis- 
tence of global weak solutions and more background discussions on the density-dependent 
Navier-Stokes equations. For the Navier-Stokes equations, it is necessary for the external 
forces to be in the functional space L 2 (Q x (0,T)) in order to obtain the global exis- 
tence. However, for the system (|l.ip - (|1.5p the term — f (u — v)pfdv does not have enough 
regularity. To overcome this difficulty, we decompose the term into two components: 

-/ (u - v)p/dv = -pu [ fdw + pl v/dv, (1.10) 

jr 3 Jr 3 Jm 3 

and we can view p J R3 v/ dv as the external force of the Navier-Stokes equations. As we 
shall see later on, the work of internal forces pu J* K3 / dv appears on the left side of Navier- 
Stokes equations and has to be nonnegative to keep the energy inequality. Meanwhile, 
we introduce a regularization function Rs (see Section 3 for the definition) as in [14] to 
construct an approximation of (|1.10|) : 

-Rs (u - v)pfdv = -pRg( / fdv)u + pR s v/dv. 

Jm. 3 Jm. 3 Jr 3 

To keep a similar energy inequality for the approximation scheme, we need to add the 
regularized acceleration in the Vlasov equation too. Then we see that the external force 
term is in L 2 (0, T;Q) and the internal forces is finite when 5 is fixed, hence we can solve 
the regularized Navier-Stokes equations. Indeed, we can obtain the smooth solution of the 
regularized Navier-Stokes equations when the initial data is good enough. The uniqueness 
and existence of the Vlasov equation can be obtained when (p,u) is smooth, see [71 HI], 
The next step is to pass the limit to recover the original system from the approximation 
scheme. We shall see that the LP regularity of velocity averages [8] and fine compactness 
of the system guarantee the existence of global weak solutions to (|l.ip - (|1.9|) . 

We organize the rest of the paper as follows. In Section 2, we deduce a prior estimates 
from (|l.l|) - (|1.5p . give the definition of weak solutions, and also state our main results. 
In Section 3, we construct an approximation scheme to (|l.l|) - (|1.5p . establish its global 
existence, use the uniform estimates and LP average velocity lemma to recover the original 
system. 



2. A Priori Estimates and Main Results 

In this section, we shall derive some fundamental a priori estimates and then state our 
main results. These estimates will play an important role in the compactness analysis 
later since they will allow us to deduce the global existence upon passing to the limit in 
the regularized approximation scheme. We shall develop these a priori estimates in the 
three-dimensional space, but they all hold in the two-dimensional space. 

First, roughly speaking, (jl.ip and the incompressibility condition mean that the density 
p(t,x) is independent of time t. In fact, we take any function j3 6 C 1 ([0, oo); fi), multiply 
(II. ip by /3'(p), use the incompressibility condition, and integration by parts over f2, then 
we have 
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Applying the maximum principle to the transport equations (jl.ip and (|1.3p . one deduces 
that 

IIpIU 00 < IIpoIU 00 , 

and also p > 0, so we have 

< p(t,x) < \\po\\l<*> (2.1) 

for almost every t. 

We now multiply (|1.2p by u and integrate over f2, and use (jl.ip . (jl.3p . and f)1.5|) to 
deduce that 



d_ 
~dt 



/ p\u\ 2 dx + 2 \Vu\ 2 dx = -2 / p/(u - v) • u dv <ix. (2.2) 
Jn Jn Jn Jr 3 



On the other hand, we multiply the Vlasov equation (|1.4p by ^j-, integrate over £1 x 
and use integration by parts to obtain 

d f f .,,■), „ f , s |v 



f w\ 

/ / /l v l dvdx- / (v • v)—pf dvdx 
at Jn JR 3 JdnxR 3 1 ^2 3) 

= -2 1 i pf\u-v\ 2 dvdx + 2 / pf(u - v)u dvdx. 
Jn Jr 3 Jn Jr 3 

We can rewrite v* as follows 

v* = Rv, where R = I — 2vv T , 

where v is the outward unit normal vector to dQ. for all x 6 d£l. By direct computation, 
one obtains that |v*| 2 < |v| 2 . Thus, we can treat the boundary term in (I2.3P as follows: 

f |v| 2 f |v| 2 f Ivl 2 

/ (v-u)—pfdvdx= (v -v)—pf dvdx + / (v • v)— pf dv dx 

JdnxR 3 z Jvu>0 1 Jvu<0 1 A \ 

i- i 1 9 ,■ i *i2 \ ZA ) 

< / (v • v)——pf dv dx — / (v* • i^) — - — pf dv* dx = 0. 

Jvu>0 2 Jv*u>0 ^ 

It is easy to get 

~T f f f dvdx — f (v • v)pf dv dx = 0, 
dt JnJR 3 JdnxR 3 

we treat the boundary term as follows 

/ (v • v)pf dv dx = / (v • v)pf dv dx + / (v-v)pfdvdx 

JdnxR 3 Jvv>0 Jvu<0 

= / (v-v)pfdvdx— / (v* • u)pf(t, x, v*) dv* dx = 0, 

Jv-v>0 Jv*u>0 



which implies that the conservation of mass: 

d_ 

It 



f dvdx = 0. (2.5) 

n Jr 3 
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This together with (|2.2p . (|2.3p and (|2.4p . we obtain the following energy equality for the 
system (fLlD - (fT3]l : 



(2.6) 



4- [ p\u\ 2 dx + ^- [ [ /(! + |v| 2 ) dwdx + 2 ( [ pf\u-v\ 2 dvdx 
at at Jci Jm. 3 Jn Jm 3 

+ 2 1 |Vu| 2 dx < 0. 
Jn 

Integrating (|2.6p with respect to t, we obtain for all t, 

[p\u\ 2 dx+[ f f(l + |v| 2 ) dvdx + 2 f f I pf\u - v| 2 dvdxdt 
Jn Jn Jm 3 Jo Jn Jm. 3 

+ 2 I f |Vu| 2 dxdt (2.7) 
Jo Jn 

|2 



< 



l m 0| 



dx+ [ [ / (l + |v| 2 ) dvdx. 
Jn Jm 3 



In Po 

By (|2.6p . it is easy to find that the global energy is non-increasing with respect to t: 
d 



Assume 



then 



p\u\ 2 dx + / /(l + |v| 2 ) d-vdx ) < 0. 
Jn Jm. 3 ) 

f\^o\_ dx+ f I f Q (i + \v\ 2 ) d-vdx <oo, 
Jn Po Jn Jm. 3 

Ill pf\u-v\ 2 dvdxdt<C, 
Jo Jn Jm 3 



and 



ll Vu llL 2 ((0,T)xn) < C, (2.8) 

sup \\p\u\ 2 \\ LHQ) <C, (2.9) 

0<t<T 

for any given T > and some generic positive constant C. Moreover, by the Poincare 
inequality we obtain 

II u IIl 2 (0,T;H1(O)) < C. (2.10) 
The maximum principle applied to (jl.4p implies that 

|l» < C||/o||l- (2.11) 



for all t G [0, T]. Moreover, fo > implies / > for almost every (£, x,v). Then, by the 
conservation of mass (|2.5p and (|2.1ip . one has the following estimate: 



L°°((0,T)xftxR 3 ) + ||/||.Z,° o (0,T;Z 1 (nxIl.. , 

<C(||/o|U °°((0,T)xfixR 3 ) 

Let w(t,x) be a smooth vector field in M 3 and let / be a solution to the following kinetic 
equation: 

a t / + vV x / + div v ((u;-v)/) = 0, 

/(0,x,v) = / (x,v), f(t,x,v) = f(t,x,v*) forxedn, vu(x)<0 
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in ft x IR 3 . DiPerna-Lions [7] obtained the existence and uniqueness of solution to (|2.13p 
when w is not smooth. Denote the moments of / by 

m k f(t,x)= f(t,x, v)|v| fe dv, 
JR 3 

M k f(t)= [ [ /(t,x,v)|v| fe dvdx, 
for any t G [0, T], x G ft, and integer k > 0. It is clear that 

M k f(t) = [ m k f(t,x)dx. 



We first recall the following lemma |14j : 

Lemma 2.1. Let w G L p (0, T; L^+^ft)) witfi 1 < p < oo and A > 1. A ssume that 

f G (L°° n x M 3 ) and m k f G L x (ft x M 3 ). T/ien, toe solution f of (gjgp s/iouW 
/wroe toe following estimates 

M k f < C ((Mfc/o) 1 /^ + (|/o| L » + 1)11^11^(0,^+^)))^ 

for all < t < T where the constant C depends only on T . 

We also recall the average compactness result for the Vlasov equation due to Di Perna- 
Lions-Meyer [8]: 

Lemma 2.2. Suppose 

-t^ + v • V x f n = div v (F n f n ) in V'{fl + M 3 x (0, oo)) 

where f n is bounded in L°°(0, oo; L 2 (ftxM 3 )) and / n (l+|v| 2 ) is bounded in L°°(0, oo; L x (ftx 
R 3 ) ; 1 ^| v | is bounded in L°°((0, oo) x ]R 3 ;L 2 (ft)). Then j K3 f n rj(v)dv is relatively compact 
in L«(0, T; L p (ft)) for 1 < q < oo, 1 < p < 2 and for r) such that mhrjs G L 1 + L°°, cj G 
[0,2). 

Remark 2.1. We shall use this lemma for the Vlasov equation to obtain the compactness 
of mo/ and mi/, which will allow us to pass the limit when e and 5 go to zero in the 
approximation . 

In this paper, we assume that 

Po > almost everywhere in po £ L°°(ft), 

mo G L 2 (ft), mo = almost everywhere on {po = 0}, J — — G L 1 (J7), (2.14) 



Po 



fo G L°°(ft x R ), m 3 /oGL 00 ((0,T) x ft). 



Definition 2.1. We say that (p, u, /) is a global weak solution to problem f)l. .8[) if 
the following conditions are satisfied: for any T > 0, 

• p>0, p G L°°([0, T] x ft), P eC([0,T];D>(n)), 1 < p < oo; 

• uei^T^oHn)); 

• p|u| 2 G L°°(0, T; L 1 ($7)); 

• f(t,x,v) > 0, for any (i,x,v) G (0,T) x ft x E 3 ; 
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• / G L°°(0, T; L°°(ft x M 3 ) nL^Hx 

• m 3 / G L°°(0,T; L x (Sl x M 3 )); 

• For any ip € C 1 ([0,T] x f2), such that div^ = 0, for almost everywhere t, 

— / mo • ip(0, x) dx + / / I — pu ■ — (pu (g) u) • Vip 
Jn Jo Jn V ^2 X5) 

+ /iVu • Vip + pp /(u — v) • ip dv dxdt = 0; 

JR3 J 

• For any <fi € C 1 ([0, T] x Q, x M 3 ) with compact support in v, such that 0(T, •, •) = 0, 

/ // / ($0 + v • + p(u - v) • V v 0) dxcivds = I [ f <f>(0,;-)dxdv; 
Jo Jn Jr 3 Jn Jr 3 

(2.16) 

• The energy inequality 

[ p\u\ 2 dx+ [ [ /(l + |v| 2 ) dvdx + 2 [ [ [ f\u- v| 2 dvdxdt 
Jn Jn Jr 3 Jo Jn Jr 3 

+ 2 f [ |Vu| 2 dxdt (2.17) 
Jo Jn 

< [ ^-dx+ [ [ (1 + |v| 2 )/ dvdx 
Jn Po Jn Jr 3 

holds for almost everywhere t G [0,T]. 

Our main result on the global weak solutions reads as follows. 

Theorem 2.1. Under the assumption (|2.14p . there exists a global weak solution (p,u,f) 
to the initial-boundary value problem (jl.ip - (|1.9p for any T > 0. 

Remark 2.2. The same existence of global weak solutions holds also in two-dimensional 
spaces. 

3. Existence of Global Weak Solutions 

In this section, we are going to prove Theorem 12.11 in two steps. First, we build a 
regularized approximation system for the original system, and solve this approximation 
system. Then, we recover the original system from the approximation scheme by passing 
to the limit of the sequence of solutions obtained in the first step. 

3.1. Construction of approximation solutions. For each e > 0, we define 

9 e := e 3 9 g 

and denote 

u e := u * 9 £ 

where 6 is the the standard mollifier satisfying 



6 G C°°(R 3 ),6 > 0, / 9dx = 1. 

Jr 3 
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By (|2.ip , all values of the solution p are bounded uniformly. The regularity of the term 
— J R 3(u — v)pfdv is not enough to solve the Navier-Stokes equation directly. Inspired by 
the work of [H], we introduce the following regularization function 

Rs(m f, mi/) = — — 1 , — F~r7' for ^ fixed 5 > °- 
l + 5J w3 fdv + S\J M3 fvdv\ 

Clearly 

< R s (m f,mif) < 1 

for any 5 > 0, and 

R s (m f,m 1 f) ->■ 1 

as 5 — >■ 0. For any fixed 5 > 0, as mentioned in the introduction, the regularized force term 



pRs / (u - v)f dv 



consists of two terms: 

p ( Rs I f dv I u and p [ Rs I vf dv 



d(pu ) 

at 



the first one is viewed as the work of internal force, and the second one is viewed as the 
external force. The regularized external force is in L 2 ((0,T) x fi), which ensures that 
the regularized Navier-Stokes equations with the work of internal force have a smooth 
solution. To keep a similar energy inequality for the approximation scheme, we need to 
regularize the acceleration term as 

Rs(u - v)pf 

in the Vlasov equation. Thus, we consider the following approximation problem: 

p t + div(pu £ ) = 0, (3.1) 

+ div(/m e ® u) - pAu + Vp + p\Rs j f dv \ u = p (r s J vf dv^j , (3.2) 

divu = 0, (3.3) 

^+v-V/ + div v (i2 5 (u-v)p/) = 0. (3.4) 

To define u e well, we need to set 

fl e = {x e ft, dist(z,dfi) > e} 

for any e > if fi is smooth. Otherwise, we can choose a smooth connected domain tt e 
such that 

{x G fi, dist(x, dfl) >e}cn e cl^cO. 

We let u e to be the truncation in Q e of u, and we extended it by to O. We define 
u £ = u £ * 0£. It is easy to find that u e is a smooth function with respect to x, and 

u e = on dtt and divu e = in R d . 

To impose the initial value for our approximate system, we need the following elementary 
variant of Hodge-de Rham decomposition (see [15]): 
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Lemma 3.1. Let N > 2, p G L°°(R Ar ) such that p > p > almost everywhere on M. N for 
some p G (0, oo). TTien i/iere exists two bounded operators Pg, Qs on L 2 (R N ) such that for 
all m € L 2 (IR Ar ), (m p ,mg) = (P p m, Q p m) is the unique solution in L 2 (R N ) of 



Furthermore, if p n G L°°(R'' V ), p < p n < p almost everywhere on M /or some < /? < 
p < oo and p n converges almost everywhere to p, then (P Pn m n , Q Pn m n ) converges weakly 
in L 2 (M> N ) to (P p m,Q p m) whenever m n converges weakly to m. 

We are ready to discuss the initial conditions for the approximation scheme (|3.ip - (|3.4p . 
Before imposing initial data, we have to point out that the initial density may be vanish 
in a domain: an initial vacuum may exist, and then in this case we cannot directly impose 
initial data on the velocity u. To remove this difficulty, we adopt the idea from [15] to 
define 



m = m p + m 



(-A)- 1 / 2 div(p- 1 m p ) = 0, (-A) 



/2 rot(m q ) = 0. 




if x is in 
if x is in Q c 



define 



Me 



Po *0e\n, 



po 2 * 0e\n, 




where d = dist(x, d£l). 
Now we define 




(3.5) 



(3.6) 



pu\ t=0 = m £ , 



and 



mg = (mopo 1/2 UpI /2 )s e c^(n). 



It is easy to see 



m e ->■ m in L 2 (Vt), m e (p £ ) 1/2 -»• m p 1/2 i: 
Relying on Lemma l3~T| we decompose as 

m^ = p £ u £ + Vq £ , where u £ , q £ G C°°(0) 
divup = in Q, 



in L 2 (n). 
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and then 

div (^( v 9o - = in Q. 



p £ Q u £ + Vq £ , where u £ G C^(Q), 



\ u o — Ug|| < e, divug = in Q. 



(3.7) 



mg->m inL 2 (0), m £ Q (p £ )- 1/2 -> m p 1/2 in L a (n). (3i 



i £ o), 



Letting 

We have 
Thus 

pu| t=0 = m £ Q = + Pq( u o - a o 
and we can impose the initial condition of u as 

u| t =o = ug. (3.9) 

Finally, we impose the initial condition for / as 

f\t=o = fo- (3.10) 

We now state and prove the following existence result. 

Theorem 3.1. With the above notations and assumptions, there exists a solution (p, u, /) 
of (I3.ip - (I3.4I) with the initial conditions (I3.5D . (I3.9P and (13.101) . and t/ie boundary condi- 
tions (dSJ), sucft i/ia£pG C°°(fix[0,oo)), u e C°°(12x[0,oo)) and / G L°°(OxM 3 x [0, oo)). 

Remark 3.1. Our approximation scheme is inspired by Lions' work on the density-dependent 
Navier-Stokes equations |15| and Hamdache's work on the Vlasov-Stokes equations |14j . 



Remark 3.2. If the initial data /o is smooth enough, we can show that the solutions are 
classical solutions. In fact, we can also show the uniqueness of such solutions. 

Proof of Theorem \3.1\ We define M as the convex set in 

C([0,T] x n) x L 2 (0,T;H^(n)) 

by 

M= {(p,u) e C([0,T] x n) x L 2 (0,T;H%(n)), 

£ < p < Co in [0, T] x f2, divu = almost everywhere on (0, T) x f2, 

ll a lli 2 (0,T;H ( 5(f7)) ^ -^}> 

where > is to be determined. Here we define a map T from M into itself as 

T(p,u) = (p,u). 
As a first step, we consider the following initial-value problem: 

^ + div(u e p) = 0, p\ t=Q = p £ , (3.11) 

in (0,T) x f2, where u e = u*# e . The construction of u e implies that u e G L 2 (0, T; C°°(Q)), 
and divu £ = in (0, T) x fi. The solution of (|3.1ip can be written in terms of characteristics: 

^- = u £ (X,s), X(x;x,t)=x, tG[0,T]. (3.12) 
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By the properties of u e G L 2 (0,T]C°°(Q)), and the basic theory of ordinary differential 
equations, we know that there exists a unique solution X of (|3.12p . Therefore, we have 

p(t,x) = p £ (X(0;t,x)), for all t G [0, T], 

It is clear that e < p < C in [0,T] x ft. Since u £ G L 2 (0, T; C°°(fi)), then p(t,z) lies 
in C([0,T];C°°(n)). By (fHTTTT) and the properties of u £ , we have §f G L 2 (0, T; C°°(fi)). 
Thus, p and ^ are bounded in these spaces uniformly in (p, u) G M. In particular, by 
the Aubin-Lions lemma, the set of p built in this way is clearly compact in C([0,T] x fi). 
The second step is to build u by solving the following problem: 

+ ^ Vu - Au + Vp + P (* I f *r) u = , (a, I v/ *,) , (3 
divu = 0, u|t=o = Uq, divup = 0, 
in(0,T) x n. Let 

e = R 5 [ fdv>0, g = R 5 [ vf dv. 



Multiplying u on both sides of (|3.13p . one obtains the following energy equality related to 

dt / -p\u\ 2 dx + / |Vu| 2 dx+ / ep\u\ 2 dx = / pgudx. 
Jn 2 Jq Jq Jq 

The right-hand side of above energy equality is bounded by 



/ pgu dx < ( p\g\ 2 dx)? • ( / p|u| 2 dx)^ 
Jn Jn Jn 



< IIPo||£ o (f7) ||5 , llL2(n)llv / P u llL 2 (n)- 
In conclusion, we obtain for all t G (0, T), 



\- j p\u\ 2 dx+ [ [ |Vu| 2 dxdt + j J ep\u\ 2 dxdt 
2 Jn Jo Jn Jo Jn 

<C I [ \g\ 2 dxdt + C [ [ p\u\ 2 dxdt + \ [ dx. 

Jo Jn 2 Jq 



lo Jn Jo Jn 1 Jn Po 

Applying Gronwall inequality, we obtain 

sup ||Hu| 2 || L i {Q) < C, 
te(o,T) 



l\/eP u llL 2 (0,T;f7) < C, 
l u lli 2 (0,T;H^(Q)) ^ 



(3.14) 



where C denotes various constant which depend only on T, f2, e, 5 and bounds on ||po||i°°(n)) 

l|Po|u| 2 || L i(Q). 

Rewriting (I3.13|) as follows 

c^ + ft-Vu-Au + Vp + au^ ^ 
divu = 0, u\t=o = Uq, divug = 0, 
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in (0, T) x Q, where 

cG L°°((0,T) x Q), b G L 2 (0,T;L°°(n)), a G L°°((0, T) x 0), 

ft G L°°((0, T) x fi), c>o>0. 
To continue our proof, we need the following lemma: 

Lemma 3.2. There exists a unique solution u of ()3. 15[) w'i/i i/ie following regularity: 

u£L 2 (0,T;H 2 (n))nC([0,T];H'(n)); Vp,| £ L 2 ((0,T) x O). (3.16) 
Proof. First, we multiply (|3.15p by ^ and use integration by parts over Q to obtain: 



du 



Of 



dx + 



dt 



1 



Vu\ 2 dx 



< / (|/i||ut| + |6||Vu||u 4 | + |a||u||u t |)dx 
Jn 



Using the Cauchy-Schwarz inequality and embedding inequality, one deduces that 



n 



du 



d 



dt 
(1 + 



dt 



lL°°(n) 



1 



dx + — / — |Vu| dx 



+ a 



2 \ 

L°°(Q) A 



(3.17) 



\Vu\ 2 dx, 



(3.18) 



where Ao is a constant from the Sobloev inequality. By the regularity of a, b and GronwalFs 
inequality, we deduce that 

u€L°°(0,T;Hl(n)), ^ G L 2 ((0, T) x n). 

We rewrite (|3.15|) as follows 

— Au + Vp = h — cu t — b ■ Vu — au, 
divu = 0, 

in f2 x (0,T), and u G Hq(Q). Let h = h — cut — b • Vu — au, and 

h G L 2 (0,T;Q), 

thus we have 

- Au + Vp = h, 
divu = 0. 

By the regularity of u and h, we conclude that p is bounded in L 2 ((0, T); if -1 (O)). We 
deduce that p is bounded in L 2 ((0,T) x fi) if we normalize p by imposing 



/ p dx = 0, almost everywhere t G (0,T). 



To normalize p, we refer the readers to |15^ [2Tj for more details. Also we conclude that 
u is bounded in L 2 (0, T; H 2 (Cl)) by the classical regularity on Stokes equation. Thus, we 
proved the regularity of (|3.16p . The existence and uniqueness of (|3. 15j) follows from the 
Lax-Milgram theorem, see for example [5]. □ 



NAVIER-STOKES-VLASOV EQUATIONS 13 

By Lemma I3.2[ there exists a unique solution to (|3. 13[) with the regularity of f|3. 16() . 
By the Aubin-Lions Lemma, u is compact in L 2 (0, T; Hq (f2)). This, with the help of 
compactness of p in M, implies that the mapping T is compact in M. 

To find the fixed point of map T by the Schauder theorem, it remains to find K > 
such that 

ll U llL 2 (0,T;H ( i(n)) ^ K - 

By (ECHD , we have 

ll U llL 2 (0,T;//i(n)) ^ ^ ) 

this K' only depends on initial data. Thus, we can choose K = K' + 1. 
Following the same argument of the proof of Lemma 13.21 we deduce that 

u G LP(0,T;W^(ft)), ^ G L"((0,T) x 0) 

for all 1 < p < oo. With such regularity of u, we can bootstrap and obtain more time 
regularity on u e and then on p and thus more regularity on u. 

In the third step, we would like to find the solutions to the following nonlinear Vlasov 
equation: 

df 

_L + vVx / + div v (i? a (u - v)p/) = 0, 

/(0,x,v) = fo(x,v), f(t,x,v) = f(t,x,v*), for x G v • v{x) < 0. 

where u, p are smooth functions obtained in step 2. The existence and uniqueness for the 
above Vlasov equation can be obtained as in [HlEj- 

Thus we have proved Theorem 13.11 □ 

Remark 3.3. The solutions (p,u,f) obtained in Theorem 13.11 satisfy the following energy 
inequality 

|(/ip|u| 2 dx+ / / i/(l + |v| 2 )^dv) 

at \Jn 1 Jn Jr 3 ' 20) 

+ / \Vu\ 2 dx+ / / R s pf(u - v) 2 dxdv < 0. 

The energy inequality will be crucial in deriving a prior estimates on the solutions (p, u, /) 
of the approximate system of equations. 

3.2. Pass to the limit as e — > 0. The objective of this section is to recover the original 
system from the approximation scheme (|3.ip - (j3.4p upon letting e goes to 0. Here and 
below, we denote by (p £ , u e , f £ ) the solution constructed in Theorem 13.11 

We take j3 G C(ft,M 3 ), use (J37TJ) and d3T3j) to find that J n f3(p £ )dx is independent of 
time t, that is, 

/ p(p e )dx= / P{p%)dx for all t G (0, oo). (3.21) 
Jn Jn 
Observing that (p £ ,u £ ,f £ ) satisfies (|3.20p . one obtains 

+ f |Vu £ | 2 dx + f [ R8p £ f(u £ - v) 2 dxdv < 0. 
Jn Jn jr3 
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Integrating it from to t, we have 



[L £ \u £ \ 2 dx+[ [ ±f(l + \ v \ 2 )dxdv 
Jn 1 Jn Jr 3 1 

+ [ [ \Vu £ \ 2 dxdt+[ [ [ R s p £ f £ \u £ -v| 2 dxdvdt (3.22) 
J o Jn Jo Jn jr 3 

< \ f pl\Ml\ 2 dx + \ I [ (l + \v\ 2 )f dxd V 

z Jn z Jn Jm 3 

for all t > 0. By (|3.22[) . one obtains the following estimates: 

\\ u6 \\L 2 (0,T;H^n)) < C, 

sup \\p £ \u £ \ 2 \\ Ll{n) <C, ( 3 - 23 ) 

0<t<T 

where C denotes a generic positive constant independent of e. 

By (|3.5p and (j3.6|) . we assume that, up to the extraction of subsequences, 

p £ -> p in C([0,T];L p (O)) for any 1 <p < oo. (3.24) 

We denote by u the weak limit of u e in L 2 (0, T; Hq (O)) due to (|3.23p . By the compact- 
ness of the embedding L p (f2) ^ Wq _1 ' 2 (0) for any p > 6/5, one deduces from (|3.24[) : 

pF^p inCQO,^;^- 1 ' 2 ^)). (3.25) 
This, together with (I3.23p . yields 

p £ u £ pu in £>'(((), T) x 0). 
Let a function 5 £ C([0, T]; L p (fl)) for any 1 < p < 00 satisfy g(0) = on and 

^+div( 5 u)=0 mV'((0,T)xn), 

then g = 0, which implies the uniqueness of the density /) when u is fixed. Thus we have 
proved that p is the solution to 
We now estimate m^f 6 : 



m f £ = f £ dv= f £ dv + f £ dv 

Jm 3 J\v\<r J\v\>r 

<C||/ e || L ^r 3 + l / |v| fc / e dv 

r J\v\>r 

for all k > 0. Taking r = (J K3 |v| fc / £ dv)fc+3 , we have 



3 3 

mot < C\\r\\L~ ( [ H k f £ dv) k+S + ( [ \v\ k f £ dv 
Letting k = 3, then 

1 /2 

IK/l^fi) < CGI/Il- + 1) Qf ^ |v| 3 / e dv) . 

Thanks to Lemma [2. 11 we conclude that m^f 6 is bounded in L°°(0, T; L 1 (J7)). This yields 

\\mof\\L°°(o,T;L2(n)) < C. (3.26) 
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Following the same argument, one deduces that 



|mi/ e || 3 < C. (3.27 



Since 

d(p £ u £ 



R s [ (u £ - v)/ £ dv is bounded in L 2 (0, T; L 3/2 (0)). 

-div(p £ u e ® u £ ) + Au e + Vp + / (u £ - v)/ £ a!v, 



(3.28 



Using the fact Rs < 1, we see that 

||pi? 5 m / £ U £ || L 2( 0iT;L 3/2( Q) ) 

< C||Po||L°°((0,T)xn)ll"T-o/ £ ||L°°(0,T;L 2 (C)) ' II u ^ II L 2 (0,T;L 6 (f2)) ) 

and 

||pi? 5 mi/ £ || iK3(0iT;L 3/2 (n)) < C||po||L-((o,T)xn)ll^i/ £ H Loo(0]T;i 3 (n)) - (3.29 
Observing 

pFRs [ (u £ - v)/ £ dv = p £ R s m f £ u £ - P R 5 m 1 f £ , 
and using (|3.28[) and (|3.29|) . we obtain that 



at 

and in particular, Vu £ is bounded in L 2 ((0, T) x O) and 

p £ i? 5 / (u £ -v)/ £ (iv is bounded in L 2 (0,T;L 3 / 2 (f])) 

while p £ u £ (g) u £ is bounded in L 2 (0, T; (O)), one obtains that 

d^" £ ) is bounded in L 2 (0,T;H~\n)). 

By Theorem 2.4 of [15] . we obtain 

y/fvf-ty/pu m L p (0,T;L r (n)) 
for 2 < p < oo and 1 < r < , and thus 

p £ u £ ^pu in L p (0,T;L r (^)) 

for the same values of p, r. 

Applying Lemma 12.21 to (13. 4j) . we obtain 

mof £ — > mo/, mi/ £ — )• mi/ for almost everywhere (i, x). (3.30 

By (|3.26j) and (13.27p . the relation (I3.30P can be strengthened to the following statement 

m / £ m / strongly in L°°(0, T; L 2 (tt)), 

mif -frai/ strongly in L°°(0,T;L 3/2 (ft)). 

By (Iffll), we have 

p £ m f £ -> pm / strongly in L°°(0,T;L *k (O)), (3.32 

and 

p £ m 1 f £ ^ p mi f strongly in L°°(0,T;L^+3(O)). (3.33 
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Thanks to (|3.32|) - (|3.33[) and the weak convergence of u £ in L 2 (0, T; Hq (fi)), one has 

Rs / (u £ — v)p £ f £ dv — > Rs / (u — v)p/ dv in the sense of distributions. (3.34) 

Jr 3 Jr 3 

The next step is to deal with the convergence of div v (i2^/9 e (u £ — v)/ £ ). Let 4>(v) £ V 
be a test function, we want to show 



lim ( f f (R 5 p £ (u £ -v)f £ )V v (j)dvdx ) 
£ ^° \Jn Jr 3 ) 

[ [ (R 5P (u - v)/)V v dvdx, 
Jn Jm 3 



(3.35) 



which can be reached by (|3.34p . 

We consider a test function ip £ Cg([0, T] x 0) such that div^ = 0, and a test function 
4> £ C 1 ([0,T] x Q x M 3 ) with compact support in v, such that 4>{T, •, •) = 0. The weak 
formulation associated with the approximation scheme (|3.ip - (|3.4|) should be 



and 



<n Jr. 3 

By ([MD-dSI), we have 



f Pq u o • ^(°> x)dx+ [ [ {-p £ u £ ■ d t tp - (pu £ ® u £ ) • Vip 
Jn Jo Jn 

- Vu £ -Vtp + tp- R 5 (u £ - v)p £ f £ dv} dxdt = 0; 
Jr 3 

[ [ [ f £ {d t (f> + v V x (p + R s (u £ - v)p £ ■ V v 4>) dxdvds 
Jo Jn Jm 3 

/o^(0, ■,•) dxdv. 



(3.36) 



(3.37) 



/9 u ■ tp dx = / m • — > I mo • (p dx as e — > 0, 
for all test functions 

All the above convergence results in this subsection allow us to recover f|2.15[) - f|2.16|) by 
passing to the limits in (|3.36p and (|3.37p as e — > 0. 

Prom (|3.22p . the solution (p £ ,u £ ,f £ ) satisfies the following: 

[ lpt\u £ \ 2 dx + [ [ \f £ {\ + |v| 2 ) dxdv 
Jn 1 Jn Jr 3 2 

+ 11 \Vu £ \ 2 dxdt + !!! Rsp e f e \u £ -v\ 2 dxdvdt 
Jo Jn Jo Jn Jr 3 



< \ [ pl\ul\ 2 dx + \ f f (1 + |v| 2 )/ dxdv. 
z Jn z Jn Jr 3 



The difficulty of passing the limit for the energy inequality is the convergence of the term 
Jo InxR 3 RsP e f £ \ u£ ~ V P dvdxdt. We follow the same way as in [HI [23] to treat the term 
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as follows 

/ / Rsp £ f £ \u £ - v\ 2 dvdxdt 

Jo W ^ (338) 

= / / (R 5 p £ f £ \u £ \ 2 - 2R s p £ f £ u £ v + R s p £ f £ \v\ 2 ) dxdvdt. 

Jo i!2xR 3 

By the embedding inequality, we have 

u £ ^u weakly in L 2 (0,T;L 6 (O)). (3.39) 
By ([3T32]) . ([S32D , we deduce that 

Rsp £ mof £ \u £ \ 2 -> R s pm f\u\ 2 weakly in L^O, T; ft) 
as e — > oo. Similarly, 

Rgp £ mif £ u £ — > Rsprriifu weakly in L 1 (0,T;r2) 

as £ — > 

Finally, let us look at the terms: 

/ / R s p £ f £ W\ 2 dvdxdt- Iff R 5 pf\v\ 2 dvdxdt 

Jtt Jr 3 Jo Jn Jr 3 



/ / m2f £ dxdt + C\\p\\l°° / / (^2/ — m2f £ )dxdt 
Jo Jn Jo Jn 



<C\ P £ -p\\ Lx 
= h + h- 

It is clear to see that I\ — )■ as e — > 0. For the term because 

f -± f weak star in L°°(0, T; x M 3 )) 

for all p G (l,oo] and m3/ £ is bounded in L°°(0, T; L 1 (r2)), then for any fixed r > 0, we 
have 

/ / f £ \v\ 2 dxdvdt = [ f X (\v\ <r)\v\ 2 f dxdvdt + 0(-) 
Jo JnxR 3 Jo JnxR 3 r 

uniformly in e where \ is the characteristic function of the ball of K 3 of radius r. Letting 
e — > 0, then r — > 00, we find 



f I f £ \v\ 2 dxdvdt ^ [ [ f\v\ 2 dxdvdt, 

Jo JnxR 3 Jo JqxR 3 



10 JtlxR 3 JO JQxl 

which means I2 — > as s — > 0. Thus, we have proved 



/ I R s p £ f £ \u £ -v\ 2 dvdxdt [ [ R s pf\u-v\ 2 dvdxdt (3.40) 

Jo JttxR 3 Jo JnxR 3 



as £ — y oo. 
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We observe that 



/ Pl\"l\ 2 dx= f — \m E -Vq £ \ 2 dx 
Jn Jn Po 

= I C^T + ^ " | ^ + Vg8) • Wo) dx (3.41) 
in V Po Po Po / 

where we used Lemma 13.11 
Using divug = 0, one obtains 

/ Po> | 2 dx+ [ ^j^dx = [ ^ dx. (3.42) 
Jn Jn Po Jn Po 

Letting e -)• 0, using (ET8|) . (13321 . (COOi (|3~42jl . and the weak convergence of (p e , u £ , f £ ), 
we obtain 



/ ip|u| 2 dx + [ [ L 
Jn z Jn Jr 3 z 

Vu| 2 dxdi + III Rspf\ u ~ v | 2 dxdvdt 
Jo Jn Jr 3 



/(l + |v| z ) dxdv 
ft 



<\ f^dx + l [ [ (1 + |v| 2 )/o dxdv. 
z Jn Po z Jn Jr 3 

So far, we have proved the following result: 

Proposition 3.1. For any T > 0, there is a weak solution (p s , u 5 , f s ) to the following 
system: 

p t + div(pu) = 0, 
d{pu) 



+ div(pu (g) u) — /iAu + Vp = pRs / (u — v)/ dv, 

Jr 3 



dt 

divu = 0, 

df 

_L + v . V / + div v (^(u - v)pf) = 0. 

with the initial data u(0, x) = uq and /(0, x,v) = /q(x,v), and boundary conditions 
u(t, x) = on <9f2, 

/(t, x, v) = /(t, x, v*) /or x G dQ, v • i/(x) < 0. 
In additional, the solution satisfies the following energy inequality: 

[ \p\u\ 2 dx+ [ [ 1/(1 + | v | 2 ) dxdv 
Jn z Jn Jr 3 z 

/ / |Vu| 2 dxdt+ III Rspf\u-v\ 2 dxdvdt 
Jo Jn Jo Jn Jr 3 



+ 



'o Jn Jo Jn 

<\l^dx + \ [ [ (1 + |v| 2 )/ dxdv. 
z Jn Po z Jn Jr 3 
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3.3. Pass the limit as 5 — > 0. The last step of showing the global weak solution is 
to pass the limit as 5 goes to zero. First, we let (p s ,f s ,u s ) be a solution constructed by 
Proposition 13. 11 It is easy to find that all estimates for (p e , f e , u e ) still hold for (//, f s , u s ), 
thus we can treat these terms as before. 

It only remains to show the convergence of the terms 

/ R s p 5 f S (u 5 -v)dv, and div(R sP \u 5 - v)). 

The next step is to deal with the convergence of div v (Rg(u s — v)p s f s ). Let <^>(v) £ V 
to be a test function, we want to show 



lim 

5-5.0 



Rsp s u 5 



pu 



fV v 4>dv ] dx 



n Jr 3 



Rsp 5 f s vV v ^ 
'n Jr 3 

pfv\7 v (fi dv dx. 



(3.43) 



To prove (|3.43[) , we introduce a new function Qs = 1 — Rs (see [H] ) , it is easy to see that 

Qs ->■ as 5 ->■ 0. 

Writing 



<5„<5 



/*V v <£dv dx 



pV 



.f V v c/>dv 



(3.44) 



Q 5 pV / / 5 Vv^dv) dx. 
On one hand, applying Lemma 12.21 to (|3.4p . we see that 

/ f S V v (j)dv ^ / fV v (fidv almost everywhere (i, x). 

il 3 JR 3 



It is easy to see 



rv v </>dv 



<C|m / 5 |. 



(3.45) 
(3.46) 

(3.47) 



This, combined with (|3.31|) . strengthens (|3.45|) as follows: 

/ f 5 V v <pdv^ [ /V v 0dv strongly in L°°(0, T; L 2 (0)). 
Jr 3 Jr 3 

By the convergence of p s , (|3.47|) and the weak convergence of u 5 in L 2 (0, T; Hq (fl)), one 
deduces 

/ 5 V v 0dv J dx- 



pV 



f2 

On the other hand, 



pu 



fV v (f> dv dx 



f Q s p s u 5 ( [ f s V v( f> dv) dx 
Jn \Jr 3 J 

<C [ Q s m f\u s \ dxdt 
Jn 

L°°(0,T;L 2 (n))ll u IIl 2 (0,T;L6(Q))IIQ<5||l2(0,T;L3(Q)) 
< C||Q<5||L 2 (0,T;L3(n)) ; 



(3.48) 
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which yields 
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Q s p 5 u s [ I f°V v cj>dv 



as 5 — > 0, 



where we used that m f 5 is bounded inL°°(0, T; L 2 (Sl)) and u 8 is bounded in L 2 (0, T; L 6 (S1)), 
and Q s -> strongly in L 2 (0, T; L 3 (0)). 

So we have proved the convergence of the first integral on the left of (|3.43p . We can 
treat similarly the convergence of the second integral of (|3.43|) . Thus, we finish the proof 

of dsasp . 

To complete the proof of Theorem 12.11 it only remains to check that (p, u, /) satisfies 
the energy inequality (|2.17p . In order to verify the energy inequality (|2.17|) . we need to 
show 

r-t 



R s p s f s \u s - v| a dxdvdt 



o Jn 



o Jn 



pf\u — v| dxdvdt 



(3.49) 



as 5 ->■ 0. 
Denote 



-<5j«5| u 5 _ v |2 dvdx ^ 



where 



E 



E 5 = E{ -2E 5 2 +E s 3 , 



p S f 5 \u 5 \ 2 dvdx = I P m f \vL°\ z dx 



n Jr 3 



P<5|„<5|2 



and 



Ei 



4 



p S f S u S v dvdx = I p°m 1 f°u° dx, 



p & f 6 \w\ 2 dvdx = I p 6 m 2 f dx. 



Q JR 3 



Write R s E d = E d -Q 5 E d , we consider the convergence of E d first. 
Since 



T 



o Jnx 

T 



p s f 5 \u s \ 2 dvdxdt 



pf\u\ 2 dvdxdt 



Jflx 



< 



o Jn 

T 



P° ~ p)m f s \u s \ 2 dxdt + 



T 



o Jn 



p(m f d -m f)\u s \ 2 dxdt 



o Jn 



pm f ( |u 5 | 2 - |u|^ ) dxdt, 



then 



f E{dt^ I [ 
Jo Jo Jn 



for all t > 0. Similarly, we obtain 



E 5 2 dt 



pf\u\ 2 dvdxdt as 5 — > 



p/uv dvdxdt as 5 — > 



JS! 



for all i > 0. 
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Finally, let us examine 

III p 5 f 5 W\ 2 dvdxdt- fll pf\v\ 2 dvdxdt 
Jo Jn Jr 3 Jo Jn Jr 3 

< \\p s — p\\l°° / / m2f S dxdt + C\\p\\l°° / / (^2/ — m,2f S )dxdt 
Jo Jn Jo Jn 

= h+h- 

It is clear that I\ — > as 5 — > 0. For the term I2, because 

f -± f weak star in L°°(0, T; L p (0 x R 3 )) 
for all p E (l,oo] and m^f s is bounded in L°°(0, T; L 1 (r2)), then for any fixed r > 0, we 



have 



/ / f\v\ 2 dxdvdt = j f x(|v| < r)|v| 2 / 5 dxdvdt + 0(- 

Jo Jf2xR 3 io JftxR 3 r 



uniformly in 5 where % is the characteristic function of the ball of IR of radius r. Letting 
5—7-0, then r — > 00, we find 

/ I f s \v\ 2 dxdvdt -> / / /|v| 2 dxdvdt, 

Jo JQxR 3 Jo JQxR 3 

which means ^2 — >• as (5 — > 0. Thus, we have proved 

I E s dt^ ((I pf\u-v\ 2 dvdxdt as 5^0. 
Jo Jo Jn Jr 3 

In order to show (|3,49p . it remains to show that 

Q s E s dt ^ as 5 0. (3.50) 



o 

By the Holder inequality, we obtain 
io Jn 



QsP °m Q f 6 \a°\ i dxdt 

< C|| ( 55||L2( 0>T;i 6(Q))||mo/ <5 ||Loo( 0>T;i 2(^)) \\u S \\ L 2( 0i t ; lG(Q))- 



Qsp m f\u 5 \ 2 dxdt ->0 as 5 ->■ 



This, together with the definition of Q5, implies that 

Jo Jf7 

for all t > 0. Following the same argument, it is easy to see 

ft 



Qsp m 1 f a u d dxdt ^ as <5 0. 



JJ) 



We write 



Qa/|v| 2 / 6 dvdxdi = / / / Qsp 5 \v\ 2 f 5 dvdxdt + Q S 
Jn Jr 3 Jo Jn J\v\<r r 



til x(H<r)Qsp 5 \v\ 2 f 5 dvdxdt + -Q 5 
Jo Jn Jr 3 r 
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uniformly in 6, where x( x ) is a characterized function. We have 

p 5 -»• p in C{[0,T];L p (n)) for any 1 < p < oo, 
and by the definition of Q s , we have 

x(|v| < r)Qsp S — >■ strongly in L p (0, T; L q (£l)) for any 1 < p, (/ < oo. 
It follows 

/ / / Qsp S \v\ 2 f S dvdxdt 
when letting <5 — 5- and r — > oo. Thus, we have proved that (|3.50p . and hence have proved 

Thanks to the convergence facts and the convexity of the energy inequality, we deduce 
(|2,17p from energy inequality in Proposition 13. 11 
The proof of Theorem 12.11 is complete. 
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